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Misleading variations in estimated rotational frequency splittings of 
solar p modes: Consequences for helio- and asteroseismology 



O 

(N 

<: 

in 



o 



> 

00 



o 

(N 



Anne-Marie Broomhall/* David Salabert,^'^ "^ William J. Chaplin/ Rafael A. Garcia,^ 
Yvonne Elsworth/ Rachel Howe^ and Savita Mathur^ 

^School of Physics and Astronomy, University of Birmingham, Edghaston, Birmingham B15 2TT 

^Instituto de Astroflsica de Canarias, E-38200 La Laguna, Tenerife, Spain 

^ Departamento do Astroflsica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain 

'^Universite de Nice Sophia-Antipolis, CNRS UMR 6202, Ohsen'citoire de la Cote dAzur, BP 4229, 06304 Nice Cedex 4, France 
^LaboratoireAIM, CEA/DSM-CNRS-Universite Paris Diderot, IRFU/SAp, Centre de Saclay, 91191 Gifsur-Yvette, France 

6 7 



'High Altitude Observatory, NCAR, PO Box 3000, Boulder, CO 80307, USA 



ABSTRACT 

The aim of this paper is to investigate whether there are any 1 1-yr or quasi-biennial solar cycle- 
related variations in solar rotational splitting frequencies of low-degree solar p modes. Although no 
1 1-yr signals were observed, variations on a shorter timescale (~ 2 yrs) were apparent. We show that 
the variations arose from complications/artifacts associated with the realization noise in the data and 
the process by which the data were analyzed. More specifically, the realization noise was observed 
to have a larger effect on the rotational splittings than accounted for by the formal uncertainties. 
When used to infer the rotation profile of the Sun these variations are not important. The outer 
regions of the solar interior can be constrained using higher-degree modes. While the variations 
in the low-/ splittings do make large differences to the inferred rotation rate of the core, the core 
rotation rate is so poorly constrained, even by low-/ modes, that the different inferred rotation 
profiles still agree within their respective Icr uncertainties. By contrast, in asteroseismology, only 
low-/ modes are visible and so higher-/ modes cannot be used to constrain the rotation profile 
of stars. Furthermore, we usually only have one data set from which to measure the observed 
low-/ splitting. In such circumstances the inferred internal rotation rate of a main sequence star 
could differ significantly from estimates of the surface rotation rate, hence leading to spurious 
conclusions. Therefore, extreme care must be taken when using only the splittings of low-/ modes 
to draw conclusions about the average internal rotation rate of a star. 

Key words: Methods: data analysis. Sun: activity. Sun: helioseismology. Sun: oscillations, stars: 
oscillations 



1 INTRODUCTION 

The frequencies of the Sun's acoustic (p-mode) oscillations vary 
throughout the solar cycle with the frequencies of the most promi- 
nent modes being at their highe st whe n solar activity is at its max- 
■ " 'imum (e.g. Woodard & Noves! T 985I: Mle . Regulo & R oca Corte j 
1989; Els worthetalJ ll990: Sal abert et alT l2004: Chaplin et alj|2007l : 

2yr) can 
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2007). Mid-term sign als (< 2yr) can also 
be observed in p-mode frequenci e s (e.g . iBroomhall et"aD l2009l : 
ISalabert et all l2009l : iFletcher et all l201ol) . Solar cycle variations 
also observed in other p-mode parameters such as pow- 



and lifetimes (e.g. Salabert. Jimenez-Reves & Tomczvlj l2003l : 



iJimenez-Reves et all |2004|) . With the advent of long, continuous, 
and high-quality ast eroseismic observations of solar-like stars (e.g. 
IChaplin et 

M l201lh_ by the Convec tion Rotation and Planetary 
Transits (CoRoT; [ Baglin et alj |2006|) space mission and Kepler 
teorucki et Zll201(il) ^ it is possible to study changes in p-mode prop- 
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erties ('Garcia et al."2010'; 'Salabert et al."201 \\h and to measure rota- 
tional splittings (Ballot et al. 2011: B eck et al..2012) . 

Solar rotation splits p modes into 2/ + 1 azimuthal orders 
(m). The difference in frequency between azimuthal orders of a 
mode can be used to infer the Sun's internal rotation profile. Vari- 
ations in the splittings of intermediate-/ modes (5 < / < 300) 
have been used to observe the torsional oscillation dHowd 120091 
and references therein). No noticeable 11-yr solar cycle changes 

in the p-mode splittings of low-/ data have been observed to date 
I — " ' ' II 1 

(e.g. Jimenez. Roca Cortes & Jimenez-Reyes 2002; Gelly et al. 2002|; 

Garcia etal. 2004, 2008; Salabert et al. 2011a). Howe et al. (2000) 
used the rotational splittings of medium-/ p modes to infer the evo- 
lution of the r otation rate a t the base of the convection zone through 
the solar cycle. iHowe et alj found a 1.3 yr periodicit y in the data, how - 
ever, the signal seemed to disappear after 2000 (H owe et"aDl2007h . 
Furthermore, the result wa s not confirmed by othe r investigations (e.g. 
iBasu & Anti'3l200lll2003h . lSalabert eX^ ( 1201 lij) observed mid- term 
variations in the rotational splittings of low-/ modes. Here we examine 
whether solar cycle-related changes, potentially including mid-term 
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variations, can be observed in the most up-to-date low-/ p-mode split- 
tings. We also discuss the consequences of observed variations in the 
rotational splitting on solar rotation profile inversions. 

Acoustic p modes are not the only oscillations that propagate 
through the interior of a star. Gravity (g) modes can propagate only 
in regions of stable stratification and so are trapped within the cen- 
tral regions of a star, beneath the convection zone. Mixed modes show 
characteristics of both p and g modes, carrying information on both 
the outer envelope and the central core. Mixed modes have been de- 
tected in stars that have evolved off the main sequence and the split- 
tings of mix ed modes have been used to infer the rotation rate of 
stellar cores ( iBeck et aLllioTil : iDeheuvels et"ai]|2012l, and references 
therein). However, here we consider only the inferred average inter- 
nal rotation rate of main sequence stars, like the Sun, that do not have 
detectible g modes and mixed modes. We must, therefore, rely on the 
splittings of low-/ p modes only and so any variations in the splittings 
could be important. 

The structure of this paper is as follows: In Section|2]we describe 
how the frequency splittings were obtained and discuss the presence 
of any periodicities. In Section[3]we compare the results obtained from 
the real data with those obtained from simulated data. In Section|4]the 
consequences of the variations in the splittings for inversions of the 
internal rotation profile of the Sun and other stars are discussed. The 
main results of the paper are summarized in Section |5] 



2 EXTRACTION OF ROTATIONAL SPLITTING 
FREQUENCIES 

We have analyzed the p-mode rotational splitting s observed by 
the Bi r mingham Solar-Osc illations Network (BiSON; lElsworth et alj 
Il995al : IChaplin et alj|l99^ during the last two solar cycles in their 
entirety i.e. from 1986 April 14 to 2010 April 7 The Global Oscilla- 
tions at Low Frequencies (GOLF; iGabriel et al.lll995l) instrument on 
board the ESA/NASA Solar and Heliospheric Observatory (SOHO) 
spacecraft has been collecting data since 1996 and so we have been 
able to analyze the velocity data (following the methods described in 
Ijimenez-Reve s et al. 2003 ; Garcia et al. 2QQ5) covering almost the en- 
tirety of solar cycle 23, i.e., from 1996 April 1 1 to 2010 April 7. 

The precision with which p-mode parameters, such as the rota- 
tional splittings, can be determined is directly related to the length 
of data set under consideration. Consequently, p-mode frequencies 
are often determined from data sets whose lengths are of the order 
of years. However, a compromise must be made here regarding the 
appropriate length of data set for study so that changes as the so- 
lar cycle evolves can be resolved. Therefore the observations made 
by GOLF and BiSON were divided into 182.5-day-long independent 
subsets. After 1996 April II, when both sets of data were available, 
we ensured that the start times of the subsets from each observational 
program were the same. 

Estimates of the mode splittings were extracted from each subset 
by fitting a modified Lorentzian model to the data using a standard 
likelihood maximization method, where the rotational splitting was a 
free parameter in the model. 

The unweighted mean mode splitting for a particular subset was 
then determined by averaging estimates in the frequency range 2000 - 
3 100 yuHz. The lower limit of this frequency range (i.e., 2000 yuHz) was 
the lowest frequency for which the signal-to-noise value allowed good 
fits to be obtained in all the subsets. Furthermore, below this frequency 
the m ode frequencies experien ce almost no sola r cycle shift (e.g. 
iLibbrecht & Woodard 1990; Els worth et aT]ll994 . Above 3100 yuHz 
it becomes difficult to accurately determine the mode splittings as the 
widths of the modes become larger than the splitting itself and so the 
mode components overlap in frequency, eventually becoming a single 



Table 1. Pearson's and Spearman's rank correlations between the BiSON and 
GOLF splittings 





/= 1 


/ = 2 


Pearson's 


0.84 


0.63 


Spearman's 


0.80 


0.63 



peak. The mean splittings were calculated for the / = I and / 
modes separately. 

Two different fitting codes have been used to extract the mode fre- 



quencies jSalabert et id]|2007l ; iFletcher et al.ll2009h . both giving sim- 
ilar results. For GOLF the mean absolute difference in the estimated 
average splittings was 0.3(T for / = I and 0.4cr for / = 2, and the differ- 
ence in the estimated splittings from the two codes was always within 
Icr. For BiSON, the mean absolute difference in the estimated split- 
tings was 0.5(T for / = 1 and 0.6it for 1 = 2. Differences in the way 
the fitting procedures dealt with the BiSON window function meant 
that before 1992, when the BiSON duty cycle was often below 50 per 
cent, the splittings obtained from the BiSON data were separated by 
up to 3(T. However, after 1992 the maximum separation of the BiSON 
splittings obtained by the two different fitting methods was I. Icr. For 
clarity, we onl y show the results of one method, which is described in 
iFIetcher et aI.L 

Many authors have examined the uncertainties associated 
with estimating splittings from Sun-as-a-star observations and 
the biases associated with the estimated splittings are now well 
known (e.g. |Eff-D arwich & Korzennik 1998; Appourchaux et all 
'2000'; Ch aplin et al]|2006l ; iGarciaet alj|2008l) . We follow the advice 



of iChaplin et al 



1 20061) in order to minimize these biases. For ex- 



ample, IChaplin et alj found that one of the main sources of bias in 
estimated splittings arises if the amplitude ratios of the m compo- 
nents inside a multiplet are not set correctly. We therefore fix these 
ratios a t previously determined and accepted values for B iSON and 
GOLF dChapIin et alJl2006l ; ISaIabert, Ballot & Garcfal201 ID . Further- 
more, the amplitude of the majority of known systematic errors be- 
come appreciable above approximately 3400 yuHz, which is above the 
upper limit on the frequency range we examine here. 



2.1 Variation with time of the rotational splitting frequencies 

Let us now consider whether the rotational splitting measured with the 
methodology outlined above show any significant variations with time. 
Fig.[T]shows the mean determined splittings as a function of time, a/(0, 
for GOLF and BiSON. There is no obvious evidence for an 11 -yr so- 
lar cycle e ffect on the splittings. This is consistent with previous re- 
sults (e.g. Jiminez^oc^Corti^^^imine^^^e^ 20021 ; iGellv et aP 



suits (e.g. I Jimene z. Roca Cortes & Jimenez- RevesI 120021 ; IGellv et alj 
l2002l ; lGarcia et alJl2004l2008l ; ISaIabert et al.ll201Iah . However, there 



are discernible mid-term (~ 2yr) varia t ions, w hich is also in agree- 
ment with the results of ISalabert et al ] teoiiah . There is significant 
correlation between the variations in the splittings observed in the two 
data sets (see Table[T), such that there is less than a 0.05 per cent prob- 
ability of the correlations occurring by chance. 

There is no correlation between the / = 1 and / = 2 mode split- 
tings for either the GOLF or the BiSON data. The splittings of the 
1=1 modes show more variation than the 1 = 2 mode splittings. 
This is not surprising and indicates that, in Sun-as-a-star data, the ro- 
tational splittings estimated from 1 = 2 modes ar e more sta ble than the 
splittings estimated using 1=1 modes (see e.g. IChaplin e t al. 2006). 
However, as we now show, this difference is accounted for by the for- 
mal uncertainties associated with the estimated splittings. We define 
di to be the temporal mean of a;(f) and dai(t) = ai(t) - dj. Fig.[2l which 
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Figure 1. Mean splittings observed in GOLF and BiSON data. The top panels show the results found using the GOLF data and the bottom panels show the BiSON 
results. The left-hand panels show the 1=1 results and the right-hand panels show the 1 = 2 results. 



Table 2. Standai'd deviation of the error normalized residuals (plotted in Fig. 

E). 



/ GOLF 


BiSON 


1 0.82 


1.14 


2 1.07 


1.00 



shows Sai(t)/(Ta (where cr„ is the formal uncertainty associated with 
a;(r)), demonstrates that the scatter in the plots is similar for / = 1 and 
2, indicating that the relative variation in the splittings is comparable. 
Table |2] gives the standard deviations of the error normalized resid- 
uals, 6ai(t)/aa. The standard deviations of the normalized residuals 
are close to unity for both GOLF and BiSON data, which is what we 
would expect to observe for a totally random process with appropriate 
error bars. 



2.2 Periodograms of the observed splittings 

To further investigate whether the observed variations are significant, 
we computed periodograms of the observed splittings. These peri- 
odograms are plotted in Fig. |3] In calculating the periodograms the 
data have been oversampled by a factor of 10. The cut-off frequency 
at the low end of the periodograms, below which no information can 
be obtained, is 0.07 yr"' and 0.042 yr"' for GOLF and BiSON re- 
spectively. Also plotted in Fie. [3] are the 1 per cent false alarm sig- 
nificance levels teh aplin et alj|2002h . which were determined using 
Monte Carlo simulations. 200,000 noise data sets were simulated, us- 
ing a normal distribution random number generator, to mimic those 
plotted in Fig.[T] The standard deviation of each point on the data set 
was taken as the Icr uncertainty associated with the rotational split- 
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Figure 2. Relative variation in the error nomialized residuals, (5a;(0/o"o). The 
top panel shows the results for GOLF and the bottom panel shows the results 
for BiSON. The red squares show the 1=1 results and the black crosses show 
the 1 = 2 results. 
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Table 3. Null probability of observing the peaks in the BiSON periodogram 
(Fig.|3) in noise. 



/ 


Frequency (yr ') 


Probability (%) 


1 


0.08 


0.190 


1 


0.58 


0.797 


1 


0.67 


0.679 


2 


0.92 


0.337 



tings plotted in Fig.[T] The simulated data sets were then used to cre- 
ate periodograms and the amplitudes observed in the simulated peri- 
odograms were used to define the 1 per cent false alarm significance 
levels. The 1 per cent significance level is somewhat arbitrary, but was 
chosen a priori so that the number of expected false detection s is less 
than unity (see e.g. . Chaplin et al.ll2002l : lBroomhall et alfcOlOh . 

There are no significant peaks (at a 1 per cent level) in the peri- 
odogram of the GOLF data indicating that there are no significant pe- 
riodicities in the GOLF rotational splittings. However, there are three 
significant peaks in the / = 1 BiSON periodogram and one significant 
periodicity in the BiSON 1 = 2 splittings. The significant peak in the 
BiSON / = 2 periodogram (at ~ 0.92 yr"') does not correspond to any 
of the significant peaks in the BiSON / = 1 periodogram. 

Table [3] contains the probabilities of those peaks in the BiSON 
periodogram where the probability of observing the peaks by change 
is less than 1 per cent. Although two of the Z = 1 peaks are only 
marginally significant the other / = 1 peak (at 0.67 yr"') and the 1 = 2 
peak are well below the significance level. The uncertainties associ- 
ated with the splittings would have to be increased by 10 per cent for 
the / = 2 modes and 17 per cent for the / = 1 modes before these peaks 
are no longer significant. 
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Figure 3. Periodograms of the mean splittings observed in GOLF (top panel) 
and BiSON (bottom panel) data. Each periodogram has been normalized so 
that the power is unity when the mean significance level is 1 per cent. This 
means that the / = 1 and 2 significance levels are now coincident. In each panel 
the red, thin line represents the I = \ results, the black, thick line represents the 
1 = 2 results and the blue dotted line represent the 1 per cent significance level. 



3 COMPARISON WITH SOLARFLAG DATA 

To help determine whether the signals observed in the BiSON data 
are really solar in origin (e.g. perhaps associated with the solar cycle) 
or simply an artifact of the data analysis procedure we examined data 
that were si mulated for the solar Fitting at Low Angular degree Group 
(solarFLAG: IChaplin et al .l2006h . The data were simulated in the time 
domain and were designed to mimic Sun-as-a-star observations. The 
simulated oscillations were stochastically excited and damped with 
lifetimes analogous to those of real solar oscillations. Five sets of 
solarFLAG data were examined. The input mode frequencies were 
the same in each simulated data set and were constant in time. The 
mode splittings were fixed at 0.4yuHz. However, the data simulated 
the stochastic nature of the mode excitation giving access to different 
mode and noise realizations. The simulated data were approximately 
9 yr in length but were split into subsets of 182.5 d for analysis. 

The splittings obtained from a typical solarFLAG data set 
(FLAGa) are plotted in Fig.|4] Almost as much variation is observed 
in the splittings obtained from the solarFLAG data as was observed 
in the real data. Table |4] shows the mean observed splittings and the 
maximum absolute deviation from the mean for both the real and sim- 
ulated data. The maximum absolute deviations of the BiSON data are 
slightly larger than all of the simulated solarFLAG data sets, while 
the GOLF data are in reasonable agreement with the solarFLAG data. 
Since the only source of variation in the solarFLAG data is the real- 
ization noise these results imply that the changes seen in the real data 
are more than likely due to realization noise. 



3.1 Impact of the BiSON window function 

The discrepancy between the maximum absolute deviations observed 
in the BiSON results and the GOLF and solarFLAG results can be ex- 
plained by the impact of the window functions of the respective data 
sets. The duty cycles of the BiSON subsets range from 22 to 88 per 
cent, with the average duty cycle being 69 per cent. Normally the duty 
cycle of the GOLF data is above 90 per cent and the mean duty cycle 
is 95 per cent. However, the lowest observed duty cycle of the GOLF 
subsets was 42 per cent, which coincides with when control of SOHO 
was temporarily lost, known as the "SOHO vacation". BiSON window 
functions were applied to the GOLF data such that the 182.5 d Bi- 
SON subset from which the window function was extracted matched 
in time the GOLF subset to which it was applied. Table |4] shows that 
this increases the maximum absolute deviations of the GOLF data. 
Furthermore, although not shown here, a marginally significant peak 
(at a 1 per cent significance level) was observed at ~ 0.9 yr"' in the 
periodogram of the 1 = 2 splittings. 

Various BiSON window functions were also imposed on the so- 
larFLAG data, with duty cycles covering the range observed here. In 
addition to the splittings observed in the 100 per cent duty cycle so- 
larFLAG data, Fig.|4]also shows the splittings obtained from one set of 
solarFLAG data when a 182.5-d-long BiSON window function, with a 
duty cycle of 88 per cent, was imposed upon each 182.5-d-long subset. 
This is just one example of the window functions that were imposed. 
The same BiSON window function was imposed on each FLAG sub- 
set so that the average duty cycle could be varied systematically. When 
some window functions were imposed the maximum absolute devia- 
tions were observed to increase to values similar to those observed 
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Figure 4. Mean splittings observed in one of tlie simulated solarFLAG data sets 
(FLAGa). The top panel shows the / = 1 results and the bottom panel shows 
the 1 = 2 results. The black crosses with the solid line represent the results 
when the solarFLAG data had a duty cycle of 100 per cent. The red squares 
with the dashed Hne represent the results when the 182.5-d-long solarFLAG 
subsets were given a BiSON window function with a duty cycle of 88 per cent. 



in the BiSON data. However, in other cases the maximum absolute 
deviation was observed to decrease. Table |5] shows that, although al- 
tering the window function changes the maximum absolute deviation, 
the obtained deviation is not dependent on the value of the duty cycle. 
For example, the / = 1 rotational splittings observed in FLAGa varied 
more when the duty cycle was 88 per cent than when the duty cycle 
was 47 per cent. 

This is consistent with the surmise that the observed variations 
are due to realization noise and that the window function could af- 
fect the variations in the rotational splittings both constructively and 
destructively via its interaction with the noise. Consider a period of 
time when the realization noise shifts the obtained splitting away from 
the true value. If the window function blanks this period of time (i.e. 
there is no data available) then the estimated rotational splitting will 
be close to the true value and the observed variation in the rotational 
splittings will decrease. However, the window function could equally 
favour times when the realization noise shifts the observed rotational 
splitting away from the true value. The variation in the rotational split- 
tings will then increase. 

This hypothesis is strengthened further by analysis of the peri- 
odograms of the solarFLAG data. Fig. |5] shows periodograms of the 
splittings obtained from the solarFLAG data sets when various Bi- 
SON window functions were imposed. We start by considering the 
top-left panel of Fig. |5] which shows periodograms of the / = 1 split- 
tings obtained from the FLAGc data when BiSON window functions 
were imposed such that the duty cycle of the data was 59 per cent, 
68 per cent and 76 per cent respectively. Even though the splittings in- 
cluded in the simulated FLAG data were constant, significant peaks (at 
a 1 per cent false alarm level) are observed in the periodograms. Sig- 



nificant peaks are observed when the duty cycle was both 59 per cent 
and 76 per cent but not 68 per cent. Similarly in the top middle panel of 
Fig.[5]significant peaks (at a 1 per cent false alarm level) are observed 
in periodograms of the / = 1 splittings obtained from the FLAGd data 
when the duty cycle is 68 per cent but not 47 per cent or 100 per cent. 
This indicates that the likelihood of observing a significant peak is not 
simply a function of the duty cycle. The bottom left and middle pan- 
els show that the same is true for the 1 = 2 splittings: The bottom left 
panel shows that significant peaks are observed in the / = 2 splittings 
obtained from the FLAGe data when the duty cycle was 39 per cent 
and 100 per cent but not 59 per cent. The bottom middle panel shows 
that significant peaks are observed in the 1 = 2 splittings obtained from 
the FLAGb data when the duty cycle was 76 per cent but not 47 per 
cent or 100 per cent. This shows that more important than the actual 
duty cycle itself is the manner in which the window function interacts 
with the realisation noise. The top and bottom right-hand panels of 
Fig. |5] show just how different the periodograms can look simply by 
changing the window function of the data. 

We note here that Tables|4]and|5]indicate that the mean / = 1 rota- 
tional splittings are consistently higher than the input value of 0.4 /jHz, 
because the rotational ly split components overlap in frequency giving 
rise to a positive bias dChaplin et al.ll200^ . 

The results of the solar-FLAG simulations show that significant 
quasi-periodic variations can be observed in the splittings even when 
no actual, underlying variations in the splittings are present. The Bi- 
SON window function can act to both enhance or reduce the signifi- 
cance of the variation. 



4 CONSEQUENCES FOR ROTATION PROFILE 
INVERSIONS 

Let us examine the consequences of these apparent variations in the 
splittings for both helioseismic rotation inversion profiles and astero- 
seismic inferences of the mean internal rotation rate. Each rotational 
splitting is a spatially weighted average of certain volumes of a star's 
internal rotation rate and various inversion techniques have been de- 
veloped to infer the Sun's internal rotation profile from the observed 
splittings. We now consider the consequences of the observed varia- 
tions in the splittings on inferences that can be made about the rotation 
profiles of the Sun and other stars. 

4.1 Inversions of the Sun's internal rotation profile 

We have seen that, although the mean BiSON 1=1 splitting is 
0.411 ± 0.007 yuHz, the maximum and minimum observed splittings 
are 0.537 ± 0.070 yuHz and 0.261 + 0.070 yuHz respectively. We have 
already shown that the observed variation in the splittings is just an 
artifact associated with the realization noise in the data and the pro- 
cess by which the data were analyzed. However, one may erroneously 
ascribe the splitting variations to a physical change in the rotation rate 
of a region of the solar interior. We have split the solar interior into 
4 regions; the solar core (0.00i?o ^ r < 0.20Rq), the radiative zone 
(0.20Rq < r < 0.70i?o), the convection zone (0.70i?o < r < 0.95Rq), 
and the near-surface shear layer (0.95i?o < r < LOOiJ©)- We have then 
determined by how much the rotation rate of that layer would need to 
change, relative to an average 2D Regularized Least Square s (RLS) in- 
version (^Schou . Christensen-Dalsgaard & ThompsorJll994i) of the so- 
lar rotation profile, to explain th e observed differences in the rota- 
tional splittings. The inversions jThompson et al.l[l996l : Ihowb et alJ 
2005) w ere made using 15 y rs of Global Oscillations Network Group 
(GONG, iHarvev et al| 19961) data. This analy sis is similar to that per- 
formed for the Sun bv Elsworth et aP (Il995bh and more details can be 
found in Appendix IaI 
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Figure 5. Periodograms of the mean splittings observed in a selection of the simulated solarFLAG data sets. The periodograms have been normalized so that the 
mean 1 per cent significance level (green dotted lines) is unity. The duty cycle of the BiSON window functions imposed on each subset are indicated in the legends. 
The top row shows 1=1 results, while the bottom row shows 1 = 2 results. The left and middle columns show that the power of peaks in the periodograms can both 
increase and decrease as the duty cycle decreases. The right-hand column shows that changes in the duty cycle can give rise to very different periodograms. 



Table 4. Mean splittings and maximum absolute deviation from the mean. 



Data set 



Mean splitting (juHz) 



Maximum absolute deviation (pHz) 





/ = 1 


/ = 2 


/ = 1 


1 = 2 


BiSON 


0.411 ±0.007 


0.406 ± 0.003 


0.150 ±0.070 


0.061 ±0.034 


GOLF 


0.412 ±0.008 


0.410 ± 0.004 


0.090 ± 0.082 


0.056 ± 0.021 


GOLF (BiSON window function) 


0.414 ±0.008 


0.412 ± 0.004 


0.125 ±0.073 


0.068 ± 0.038 


FLAGa 


0.424 ± 0.010 


0.396 ± 0.005 


0.095 ± 0.054 


0.048 ±0.021 


FLAGb 


0.423 ±0.010 


0.397 ± 0.005 


0.064 ± 0.053 


0.059 ± 0.019 


FLAGc 


0.426 ± 0.010 


0.406 ± 0.005 


0.096 ± 0.076 


0.035 ± 0.024 


FLAGd 


0.421 ±0.011 


0.401 ± 0.005 


0.097 ± 0.077 


0.042 ± 0.045 


FLAGe 


0.413 ±0.010 


0.408 ± 0.005 


0.069 ± 0.068 


0.054 ± 0.022 



Table 5. Mean splittings and maximum absolute deviation from the mean estimated from the FLAGa data set when different window functions were imposed on the 
data. 



Duty cycle (%) 



Mean splitting (^Hz) 



Maximum absolute deviation (//Hz) 





; = 1 


/ = 2 


/ = 1 


1 = 2 


100 


0.424 ±0.010 


0.396 ± 0.005 


0.095 ± 0.054 


0.048 ±0.021 


88 


0.426 ±0.010 


0.400 ± 0.005 


0.117 ±0.063 


0.051 ±0.021 


76 


0.416 ±0.010 


0.396 ± 0.005 


0.110 ±0.054 


0.043 ± 0.023 


68 


0.435 ± 0.010 


0.392 ± 0.005 


0.121 ±0.047 


0.055 ± 0.024 


59 


0.417 ± 0.009 


0.382 ± 0.005 


0.103 ±0.062 


0.046 ± 0.020 


47 


0.421 ± 0.008 


0.399 ± 0.005 


0.095 ± 0.056 


0.088 ± 0.023 


38 


0.414 ±0.010 


0.379 ± 0.006 


0.130 ±0.063 


0.053 ± 0.035 


22 


0.406 ± 0.012 


0.403 ± 0.008 


0.162 ±0.078 


0.138 ±0.029 
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Table 6. Variation in £1/2^ required to explain a 0. 15 p^Rz change in the rota- may, therefore, erroneously conclude that the average rotation rate of 

tional splitting of / = 1 and 2, « = 17 modes. the stellar interior is dilferent to the surface rotation rate. 



Shell 


Inner radius of 


/= 1 


/ = 2 




shell (r/Ra) 


(pUz) 


(pHz) 


core 


0.00 


2.36 


3.27 


radiative zone 


0.20 


0.50 


0.48 


convection zone 


0.70 


0.47 


0.47 


surface shear layer 


0.95 


0.49 


0.49 



5 DISCUSSION 

As with pre vious stu dies, (e.g . 

JimeneZj^BcaCort^s & Jimen ez-Revej 20021 : iGellv et al] |2002| : 
Garcia et al.l2004l2008l ;iSalabert et al.l201 lah . we find no 1 1-yr solar 



Using the above mentioned maximum and minimum measured 
splittings we find that the maximum departure from the mean BiSON 
splitting was 0.150 ± 0.070 A(Hz. Table [6] shows the change in rota- 
tion rate, d(n/27r), required to change the observed splitting of / = 1 
and 2, n = 17 modes by 0.15yuHz (results for other modes can be 
found in Appendix [At. The results indicate that the changes required 
in the rotation rate of near-surface regions to explain a difference in the 
splitting of 0.15;uHz are not unreasonably large. However, observa- 
tions of higher-/ modes rule out variations in the near-surface regions 
even of this magnitude (see Appendix iBt. Information from higher-/ 
modes dominate and constrain inversions of the rotation in outer re- 
gions of the solar interior at better than a 0.5 per cent level. Although 
variations with time in the higher-/ splittings are detectable (due to 
the torsional oscillation) th ese variations are less than 0.001 yuHz (e.g. 
iHowe. Komm & Hili[T999l) and so they are approximately 100 times 
smaller than the deviations observed here. Large changes in the core 
rotation rate are required to explain the observed change in the rota- 
tional splittings. Although the required changes are so large as to be 
unlikeljlJ such a large change cannot be ruled out with current inver- 
sions of the solar rotation profile (see Appendix |B). Even low-/ modes 
carry so little information about the conditions of the core that the 
uncertainties associated with any inversions are very large. 



4.2 The rotation profile of a solar-type star 

Asteroseismic data often consist of just a few months of observations 
in which only the low-/ modes are visible. Thus asteroseismic data 
are often similar to the 182.5 d Sun-as-a-star subsets analyzed here. 
We note here that this analysis is relevant for main sequence stars 
for which no mixed modes or g modes can be detected and used to 
aid inferences of the internal rotation rate. Therefore it is difficult to 
make inversions of the internal rotation profile of a star. However, 
lightcurves of a star, such as those observed by CoRoT and Kepler, 
can show rotational modulation due to the presence of magnetic stellar 
activity, such as spots and active regions. The modulation of the light 
curve can, therefore, be used to determine the surface rotation rate of 
a star (e.g. iMosser et al.]|2005L 120091: iMathur et al.]|2010t IB allot et alj 
I2OI ih . which can then be compared to the mean internal rotation rate 
determined from the oscillations. The sidereal rotation rate of the Sun, 
as estimated from the rotation of sunspots at the surface, is approxi- 
mately 27 d. The mean rotation rate of the solar interior, determined 
from low-1 rotational frequency splittings, is 28 d. However, if we 
were to use the extreme values of the splittings obtained from a single 
182.5 d subset, similar to those available in asteroseismic studies, we 
would find the mean rotation rate to be either 22 ± 3 d or 44 ± 12 d. One 



* The rotati on rate in the core is not well constrained through p modes (e.g 
lHowdl2009l and references therein)and so, in this paper, we assumed a rate 
similar to the rest of the radiative zone. However, we note that g modes 
would provide tighter constraints on the core rotation rate jGarcia et al.l2007t 
[Mathur et al. 2008). 



cycle variation in the rotational splittings. Discernible quasi-periodic 
mid-term (~ 2 yr) variations are present in the splittings. However, we 
show that these variations are due to realization noise. 

Many authors have discussed the bias es associated with estimat 



f'ittings from Sun-as-a-star data (e.g. Eff-Darwich & Korzennik 
lApp ourchaux et al. 2000: Chaplin et al. 2006; Garci a et af 
We have shown that, despite taking the advice of lChaolin et al 
, the realization noise has a larger effect on rotational split- 
tings than accounted for by formal uncertainties, implying that the 
uncertainties on the fitted splittings are underestimated. The elfect 
of the realization noise produces an erroneous apparent mid-term 
(~ 2 yr) signal in the observed splittings that could, potentially, look 
like the quasi-biennial solar-cycle related signal that is observed in 
the mode frequenc ies ( iBroomhall et al.l l2009l: ISalabert etaTI l2009l : 
iFletcher et al.ll2O10l) . Although we have shown that the mid-term sig- 
nal is an artifact associated with the realization noise and the data anal- 
ysis process the periodicities in the splittings may be mistakenly inter- 
preted as a physical variation in the rotation rate. Such variations are 
unlikely because they would require either very large changes in the 
core rotation or changes in the near surface that, although smaller, are 
ruled out through observations of higher-/ modes. 

It is not surprising that the splittings observed by BiSON and 
GOLF are highly correlated, since they both are observing the same 
star. In fact we expect the splittings to be correlated whether these vari- 
ations are of solar o rigin or whether, as we have shown, they are due 
to realization noise. Ijimenez-Reves et al] ( |2004|) found that there was 
a strong degree of correlation between numerous mode parameters, 
such as mode linewidths and power densities, observed in GOLF and 
BiSON data. They concluded that this showed that the estimated mode 
parameters were dominated by the same mode realiza tion noise (the 
signature of stochastic excitation). iHowe et al.l ( I2OO6I) observed that, 
after the solar-activity dependence had been removed from the fre- 
quencies, significant fluctuations in the frequen cies were hi ghly corre- 
lated between BiSON, GONG, and MDI data. lHowe et aL explained 
these correlations in terms of the stochastic nature of the mode ex- 
citation, which could be interpreted as a form of realization noise. 
The high level of correlation observed here between the GOLF and 
BiSON splittings is also due to the mode realization noise. This sup- 
ports the conclusion that the mode realization noise is responsible for 
the observed variation in the splittings, which again demonstrates how 
careful one needs to be, particularly when data have correlated noise 
sources. 

Although the spurious variations do not have a significant effect 
on inferred rotation profiles of the Sun they could be important for 
any asteroseismic inferences of the average internal rotation rate of 
main sequence stars for which g modes and mixed modes have not 
been detected. Firstly because, unlike for the Sun, we must rely on 
the splittings estimated from low-/ modes only. Secondly, because we 
often have only one relatively short data set, of the order of a few 
months, from which to extract the splittings. An additional consid- 
eration is that the accuracy of the determined splitting is dependent 
on the size of the splittings relative to th e width of the modes i.e. 
whether the modes ove rlap in frequency dBallot, Garcia & Lambert! 
I2OO6I : iBallot et al1l2008l) . We have shown that care must be taken in 
using average splittings (extracted from data sets of a few months du- 



© 0000 RAS, MNRAS 000, 000^00 



8 A. -M. Broomhall et al. 



ration) to infer the mean internal rotation rate (which, for example, 
might be compared in asteroseismic studies to the surface rotation 
rate inferred directly from rotational modulation of the light curve, 
e.g.lBallot et al. 2011). Artificial data, such as those available through 
AsteroFLAG (Chaplin et al. 2008), should be used to test for poten- 
tial biases in the results. Care should be taken when interpreting the 
observed splittings of low-Z modes in both astero- and helioseismic 
data. 



APPENDIX A: INVERSION PROCESS: FORWARD 
CALCULATIONS 

A rotation profile of the Sun was simulated, based on the average of 
2D inversion p rofiles computed over 15 yrs using a Regularized L east 
Squares (RLS, ISchou. Christensen-Dalsgaard & ThomDsoiJll994) in- 
version method and Global Oscillations Network Group data (G ONG, 
iHarvev et al.) 1 19961 : iThompson etaP 1 19961 : iHowe et alj l2005h . The 
GONG instrument make spatially resolved observations of the Sun 
and is thus able to observe much higher-/ modes than those exam- 
ined here. The GONG observations were used to build a latitudi- 
nally dependent rotation profile of the Sun between 0.5 and 1.0 so- 
lar radii (Rq). Below 0.5 Rq the rotation was taken to be constant and 
equal to the mean inferred rotation rate across all latitudes at 0.5 Rq. 
We varied the rotation rate in four different regions of the solar in- 
terior; the solar core (O.OO^o < r < 0.20Rq), the radiative zone 
(0.20/?G < r < OJORe), the convection zone (0.70^o < r < 0.95Ro), 
and the near-surface shear layer (0.95/?o < r < I.OORq). The splittings 
were calculated when the rotation rate in each region was varied be- 
tween 330 < Q./2n < 530 nHz. The 1=1 and 2 spl ittings, 6v, were 
then determined (see for example IChaplin et al1ll999l) . 

To first order, the magnitudes of the splittings vary linearly with 
the rotation rate in a given region. That linear factor depends on the 
position of the specific region in the solar interior and so can be used 
to infer the change in rotation rate in a specific region that is required 
to produce a change in splitting of 0. 15 //Hz. These values are plotted 
in Fig. lAl l and are listed in Table|6]for / = 1, n = 17 and I = 2,n = 17 
(corresponding to modes at 2561 and 2621 yuHz respectively). The re- 
sults in Fig. lAll and Table |6] indicate by how much the rotation rate 
in a particular region would have to deviate from the values obtained 
from the 15-yr 2D RLS inversions to explain the observed results. The 
required deviation in rotation rate change is similar for the radiative 
zone, the convection zone and the surface shear layer and is approx- 
imately constant across the range of n examined here. However, the 
change in rotation rate required in the core is substantially larger and 
depends on n. 

As expected from the rotation kernels, the observed splittings are 
much more sensitive to surface regions than to the core. The changes 
in Table [6] can be compared to the mean rotation rate at the surface, 
which is approximately Q./2n = 430 nHz. The rotation rate in the core 
would have to be approximately 5 times faster than the mean solar sur- 
face rotation rate to account for an increase in the splitting of 0. 15 yuHz, 
while the rotation rate in the other regions of the solar interior would 
have to change by more than 100 per cent of the mean solar value. In 
other words, even if the shell responsible for the change in splitting 
is close to the surface, a large, but not inconceivable, change in rota- 
tion rate is required to explain a difference in the observed splitting of 
0.15aiHz. 
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Figure Al. Change in rotation rate required to account for a 0. 15 yuHz change in 
the observed splitting as a function of n, for / = 1 (top panel) and 1 = 2 (bottom 
panel) modes. The different symbols coiTespond to the results obtained when 
the rotation rate in a specific region was altered (see legend). 



APPENDIX B: VARIATIONS IN THE INFERRED SOLAR 
ROTATION PROFILE 

We have combined the low-degree splittings observed in the BI- 
SON and GOLF data with a set of intermediate-degree splittings 
made from the average o f three 72 d Michelson Doppler Imager 
(MDI. IScherrer et al.l[l995i^ sets observed between 1996 and 1997. 
Each set of splittings was then inverted using 2 D RLS methods 
dSchou. Christensen-Dalsgaard & Thompson! Il994l) . Fig. IB II shows 
residuals in the inferred rotation profile at the equator once a mean 
profile has been subtracted. In regions where no information is avail- 
able (i.e. the core) the 2D RLS inversion method extrapolates a solu- 
tion. In these regions the averaging kernels, the functions that describe 
how the estimate of the solution corresponds to a spatial average of the 
underlying true solution, are not well localized. The rotation profiles 
plotted in Fig. lBll were generated using two sets of splittings measured 
by BISON and GOLF. These were chosen because they correspond to 
epochs when the observed splitting deviates noticeably from the mean. 
In each case the rotation profiles differ from the mean, and from each 
other, in the deep interior but these differences are less than the formal 
uncertainties for the inversions. Therefore, the spurious changes in the 
observed splittings do not have a significant effect on the inferred ro- 
tation profile of the Sun. This is not entirely unexpected because even 
though the (model-dependent) spatial weighting functions (kernels) of 
low-l modes do probe the solar core, they are still more sensitive to the 
surface regions. Indeed, the sound speed in the deep interior is much 
larger than near the surface, so the time the modes dwell in the core is 
small relative to the outer regions. 



© 0000 RAS, MNRAS 000, 000-000 



Variations in rotational frequency splittings 9 



0.10 



N 

X 



0.05 



CS 



0.00 



"1 — I — I — I — I — I — I — I — I — I — I — I — p 



J I L 



J I I I I I I I I I I I L 



J I L 



0.0 0.2 0.4 0.6 0.8 1.0 

r/R 

— I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



0.05 



X 0.00 - 



3 -0.05\- ■■ / 
/ 



-0.10 7 



J I I I I I I I I I I I L 



0.0 0.2 0.4 0.6 0.8 1.0 
r/R 

Figure Bl. Residuals in the inferred profile at the equator once a mean profile 
was subtracted. The top panel shows the results obtained using the splittings 
observed in a 182.5 d subset that began on 2001 April 10 while the bottom panel 
shows results using the spHttings observed in the 1 82.5 d subset that began 
on 2008 October 7. The black solid line shows the results obtained using the 
BiSON data, while the Icr uncertainties are denoted by the black dotted lines. 
The red thick line shows the results obtained using the GOLF data, while the 
Icr uncertainties are denoted by the red dashed lines. 
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